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Abstract 

We study the hydrogen storage properties of planar boron sheets and compare them to those 
of graphene. The binding of molecular hydrogen to the boron sheet (0.05 eV) is stronger than 
that to graphene. We find that dispersion of alkali metal (AM = Li, Na, and K) atoms onto 
the boron sheet markedly increases hydrogen binding energies and storage capacities. The 
unique structure of the boron sheet presents a template for creating a stable lattice of strongly 
bonded metal atoms with a large nearest neighbor distance. In contrast, AM atoms dispersed 
on graphene tend to cluster to form a bulk metal. In particular the boron-Li system is found to 
be a good candidate for hydrogen storage purposes. In the fully loaded case this compound can 
contain up to 10.7 wt % molecular hydrogen with an average binding energy of 0.15 eV/H2. 
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Introduction 

Hydrogen is an abundant, clean, and renewable energy carrier.^ An important barrier prevent- 
ing the large scale use of hydrogen is storing it densely and safely under moderate conditions.^ 
Also, fast loading and unloading of the storage system is still a challenge. Storage of hydrogen 
in molecular form may be more beneficial than storage of atomic hydrogen in chemical hydrides. 
Complex metal hydrides, for instance, are often either too stable or too unstable^ and require sub- 
stantial doping to tune their stability. 5 6 Moreover, the formation and decomposition of chemical 
hydrides typically involve complicated solid-state chemical reactions, which hamper the kinetics 
of hydrogen loading and unloading. In contrast, materials that physisorb molecular hydrogen, 
such as graphene, carbon nanotubes and fullerenes, clathrates, zeolites, and metal organic frame- 
works (MOFs), are capable of achieving fast hydrogen kinetics. The binding energy of hydrogen 
to these host materials is however small, which leads to unfavorable operating conditions of very 
low temperatures or high hydrogen pressures. 7 

Recent studies show that metal doping enhances the strength of binding between hydrogen 
molecules and physisorption materials.^ An ideal metal dopant should strongly bind to the host 
material, and it should bind hydrogen molecules effectively. For instance, for storage of molecular 
hydrogen at 30 bar under ambient temperature, and delivery at 1.5 bar, the optimum adsorption 
enthalpy needs to be 0.15 eV/H2. 7 Upon deposition, the geometrical positioning of metal atoms 
over the substrate is also important, since the number of H2 molecules per dopant site should be 
maximized. To prevent steric hindrance, this usually means that the dopant atoms need to be far 
apart. Therefore, the material should resist the clustering of metal atoms. At least the binding 
energy of a metal atom to the substrate should then be larger than the cohesive energy of the 
bulk metal. This is extremely difficult to attain with transition metal dopants, since the cohesive 
energies of bulk transition metals are high. Alkali metals have much lower cohesive energies and 
are therefore more suitable as dopants. 

Recent studies addressed the alkali metal doping of physisorption materials such as graphene,! 10 ! 11 ! 12 ! 
fullerene, ! 13 ! 14 ! 15 ! carbon nanotubes ; l 1 6 1 1 7 1 1 8 I S I and MOFs 1 l 22 l 23 l Hydrogen molecules are trapped 
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around the dispersed doping metal atoms through electrostatic and polarization interactions.^ In 
addition to increasing the binding strengths to hydrogen, alkali metal doping also increases the ef- 
fective surface area available for absorption in some cases. Hydrogen adsorption then goes beyond 
a single monolayer, increasing the volumetric and gravimetric hydrogen densities accordingly. 

In order to obtain a sufficiently high gravimetric density, not only the alkali metals need to 
be lightweight, but also the substrates need to be lightweight materials. In the past decade, car- 
bon nanostructures of various dimensions have been studied extensively both theoretically and 
experimentally for various applications, including as substrates for hydrogen storage. The recent 
realization of the structurally most simple carbon network, graphene, has attracted a great deal of 
interest.^^ The lightest element that can form extensive covalently bonded structures is in fact 
boron. Boron, having one electron less than carbon, forms icosahedral clusters ,1221281 an( j j ias \ on g 
been thought not to form two-dimensional (2-D) structures similar to graphene. 

However, it has very recently been shown that the stability of planar 2-D hexagonal boron 
sheets is enhanced markedly by inserting additional boron atoms in the centers of hexagons of the 
honeycomb structure; see Figure 1 l 29 l 3 ° l 31 l 32 l These additional boron atoms modify the electronic 



structure of the boron lattice. 2-D hexagonal boron is an electron deficient system with part of the 
sp 2 bonding states unoccupied. The electrons provided by the additional boron atoms occupy these 
states, increasing the stability. 29 These new 2-D boron structures also provide insight into the un- 
usual stability of other boron-only structures with different dimensions, such as fullerenesJ^EMl 
nanotubes ; l 32 l 36 l 37 l and nanoribbons.^ 

Here, we investigate the hydrogen storage properties of the boron sheet structure that is iden- 
tified as the most stable form, the so-called a-structure. 29 We first study the interaction of H2 
molecules with this boron sheet. Next, the dispersion of alkali metal (AM = Li, Na, and K) atoms 
onto the boron sheet is considered. We study the binding and the mobility of these AM atoms on 
the boron sheet and the hydrogenation energies of these systems. All of the three AMs are found 
to bind strongly to the boron sheet. Moreover, the unique structure of the boron sheet allows for a 
stable lattice of AM atoms with a large nearest neighbor AM— AM distance, which consequently 
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Figure 1: The most stable form of the boron sheet. The primitive cell is indicated by dotted 
red lines. It contains eight boron atoms, shown as purple spheres. The irreducible part of the 
Wigner— Seitz cell (WSC) is indicated by a triangle. 

creates space for hydrogenation. The AM atoms accumulate a net positive charge, which not only 
stabilizes this structure but also strengthens the binding to hydrogen molecules. All three AMs 
lead to a higher hydrogen binding energy, but Li is by far the best doping element. 



Computational Details 

We perform first-principles calculations based on density functional theory (DFT), using the Vi- 
enna ab initio simulation package (VASP).^^ The generalized gradient approximation (GGA) 
in the form of the PW91 functional is used to approximate exchange and correlation. 41 Although 
dispersion (van der Waals) interactions are not captured by this functional, it should give a good 
description of chemical bonding and of electrostatic interactions, which are the most prominent 
interactions in alkali doped systems; see also ref.^ Binding energies of H2 molecules to similar 
systems obtained with the PW91 functional are within ~20 meV/H2 of those obtained with the 
PBE functional.^ Therefore, we only give the PW91 results in the following. 

The basis set is constructed according to the projector augmented wave (PAW) method. E2143J 
The following electrons are treated as valence: H Is, B 2s2p, Li ls2s, Na 3s2p, and K 4s3p. A 
cutoff energy of 400 eV is used for the plane-wave basis. Brillouin zones (BZs) of all structures 
are integrated with fc-point spacings of ~0.01 A _1 JSlFor instance, the BZ corresponding to the 
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(b) 



Figure 2: The most stable geometry of H2 adsorbed on a boron sheet shown in (a) top, and (b) 
side views. The molecule is positioned over the center of an empty boron hexagon with its axis 
orthogonal to the boron plane. 



unit cell shown in Figure Figure 1 is then integrated using a 12 x 12 regular grid. Gaussian 
smearing is used with a smearing width of 0.01 eV. Periodic images of the sheet along the surface 
normal are well separated by a distance of ~20 A, so that the interactions between these images 
are negligible. Atomic and geometric relaxations are carried out by employing the conjugated- 
gradient (CG) algorithm. All atomic positions and cell parameters are relaxed. The relaxations 
are assumed to be complete when the total force remaining on each atom and the stress tensor 
components acting on the cells are less than 0.02 eV/A and 2.5 kbar, respectively. Total energies 
then are converged to within 1 meV/unit cell at least. The tetrahedron method is then used to 
calculate the total energies of the relaxed structures.^ 



Results and Discussion 



Hydrogen Adsorption on Boron Sheets 



The most stable structure of a boron sheet is shown in Figure 1 It is a graphene-like structure, 
where two-thirds of the hexagon centers are occupied by additional boron atoms. The primitive 
cell of the sheet then contains eight boron atoms. It leads to a unique planar structure with a lattice 
of empty hexagons spaced at a distance of 5.06 A of one another. 

To model the interaction of hydrogen molecules with the boron sheet, we use a 2 x 2 supercell. 
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The distance between hydrogen molecules in neighboring cells is then 10.12 A, which is suffi- 
ciently large to neglect the interaction between these molecules. We consider a number of possible 
positions of the H2 molecule on the boron sheet. Since the latter contains a lot of symmetry, it 
is only necessary to consider the irreducible part of the Wigner— Seitz cell (WSC), which is indi- 



cated in |Figure Tj We scan the edges of the irreducible WSC by placing a H2 molecule at several 
positions, which are separated by 0.8 A, as measured from the H— H bond centers. At each of 
these points, we consider three different molecular orientations. In two of these orientations, the 
H2 molecular axis is parallel to the boron plane. The two orientations differ by an in-plane angle 
of 30°. In the third configuration, the H2 molecular axis is perpendicular to the boron plane. The 
H— H bond centers are placed at a distance of at least 2 A from the boron sheet. 

In total, we have constructed a set of 36 distinct initial geometries. All of the structures in this 
set are then subjected to relaxation without any symmetry constraints. These calculations reveal 
that H2 molecules tend to occupy positions over the open boron hexagons. The minimum energy 
configuration is the one in which the molecule is over the center of an open hexagon with its axis 



perpendicular to the boron plane (Figure 2 ). The distance between the boron plane and the center 
of the molecule is 3.27 A, and the calculated value of the binding energy is 47 meV. The energy 
differences between different orientations of the H2 molecule are generally small. A H2 molecule 
at this site with its axis parallel to the boron plane is only ~ 5 meV higher in energy. The potential 
energy surface obtained by varying the position of the H2 molecule over the plane has a corrugation 
of 61 meV. 

Interestingly, the interaction of H2 with the boron sheet shows similarities to its interaction 
with graphene. We consider a graphene supercell similar to that of boron, which has of course the 
typical graphene structure without atoms in the centers of the hexagons. The calculated binding 
energy of H2 to graphene is then 25 meV. Likewise, the geometry of the adsorbed H2 molecule is 
very similar to that on the boron sheet, and the distance to the graphene plane is 3.24 A. Our results 
on the graphene— H2 system are in accordance with earlier studies J ^ 47 * 48 ! We note that, although 
the binding of H2 to the boron sheet is slightly stronger than that to graphene, it is still too weak 
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Figure 3: Optimized structure of a boron sheet doped with Li atoms (yellow spheres), shown 
from (a) the top and (b) the side. Boron sheets doped with other lightweight AMs have a similar 
geometry. 

for a practical hydrogen storage application. 



Doping Boron Sheets with Alkali Metals 

Next, we consider deposition of AM atoms onto the boron sheet, with AM = Li, Na, or K. To have 
a stable material at a reasonable temperature, it is essential that the metals strongly bind to the 
boron template. In addition, meaningful hydrogen gravimetric densities can be reached only if the 
AM atoms are separated far enough from each other so that each atom has sufficient space to bind 



hydrogen molecules. Placing a single AM atom in the unit cell shown in Figure 1 leads to a lattice 
of AM atoms with a cell parameter of 5.06 A. We vary the position of the AM atom on the boron 
sheet and relax all structures, as discussed in the previous section. The optimal position of all three 



lightweight AMs turns out to be over the center of the empty boron hexagon, as in Figure 3 Upon 



relaxing, the unit cell parameters of the doped systems do not change noticeably. The optimized 
cell parameters are 5.08, 5.07, and 5.07 A for the Li, Na, and K doped boron sheet, respectively. 
The binding energies of the AM atoms to the boron sheet are calculated according to 



AMbulk "I" ^boron ^boronAM 



-E h 



(1) 
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Table 1: Calculated Binding Energies of Alkali Metals on a Boron Sheet ( |Figure 3 ) with Respect 



to Bulk Metal (i+X Isolated Atoms (E at ), Equilibrium Distances to the Boron Plane (d), and Bader 



AM 


£ b (eV) 


£at(eV) 


d(A) 


Q(e) 


Li 


0.27 


1.89 


1.54 


+0.85 


Na 


0.02 


1.11 


2.26 


+0.60 


K 


0.34 


1.21 


2.75 


+0.47 



where ^AMbuik is the total energy per atom of the AM bulk metal, £boronAM is the total energy per 
AM atom of the AM doped boron sheet, and i+oron is the corresponding total energy of the boron 



sheet. The results are shown in Table 1 For comparison, the binding energies E at , calculated from 
eq[T]by replacing ^AMbuik with the energy of an isolated AM atom, are also given in this table. The 
numbers show that all three systems are thermodynamically stable with respect to separation into 
bulk metal and a boron sheet. 



The equilibrium distance between the AM atoms and the boron plane, listed in Table 1 shows 
an increase going from Li to K, as might be expected. Electrostatic interactions are dominant in 
binding hydrogen molecules to the AM atoms The Bader charges decrease going from Li to 
K, which at first sight may seem surprising as the electronegativity of the AMs decreases along 
this series. However, comparing the nearest neighbor distances between the metal atoms on the 
boron surface with those in the bulk metals, one observes that the former are larger by 70, 40, and 
10 % for Li, Na, and K, respectively. This means that the Li atoms on the boron sheet are relatively 
isolated, but the K atoms are still in relatively close contact. In this configuration, a Li atom is 
bonded to the boron sheet and the binding has a large ionic contribution, whereas a K atom is 
bonded to the boron sheet and to other K atoms with a mix of ionic and metallic bonding. Clearly, 
the binding of a Na atom is intermediate between these two situations. 

In principle, one might increase the distance between the K atoms by decreasing their concen- 
tration on the boron sheet. However, K is already the heaviest element of the AM series considered 
here, and a reduction of their concentration increases the number of boron atoms per metal atom. 
This would seriously decrease the gravimetric hydrogen density, since the amount of hydrogen 
adsorbed depends on the metal atoms only. In conclusion, since Li is the lightest element, binds 
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Table 2: Calculated Binding Energies with Respect to Bulk Metal (E^) and Bader Charges (Q) of Li 
on a Boron Sheet at Various Concentrations, Represented by the Nearest Neighbor Metal— Metal 
Distance d nn . 



J nn (A) 


£b(eV) 


Q(e) 


10.16 


0.45 


+0.88 


5.08 


0.27 


+0.85 


2.93 


0.16 


+0.75 



strongly to the boron sheet, and discharges almost completely upon deposition, it is the best can- 
didate as a dopant for hydrogen storage purposes. 

Since the Li atoms on the boron carry an effective positive charge, one may expect a net re- 
pulsion between these atoms. We consider the binding energy E^, eq[TJ and Bader charge Q as a 



function of the Li concentration. The results are given in Table 2 In the concentration we have 
considered so far, the Li— Li nearest neighbor distance is 5.08 A. Doubling the boron surface unit 
cell in each direction and occupying only one empty hexagon by a Li atom decreases the concentra- 
tion by a factor of 4 and decreases the Li— Li nearest neighbor distance by a factor of 2. Occupying 
not only all empty, but also all filled boron hexagons increases the Li concentration by a factor of 
3. The Li— Li nearest neighbor distance is then multiplied 1 / v(3). The binding energy indeed 
decreases with increasing concentration, confirming the net repulsion between the Li atoms. There 
is a corresponding slight decrease of the Bader charge on the Li atoms. Note however that, at all 
concentrations studied, the system is stable against phase separation into a sheet carrying a lower 
concentration of Li and bulk Li. The Li concentration can therefore be controlled by the amount 
of Li dispersed on the boron sheet. 

It is interesting to compare the Li doping of the boron sheet to that of a similar system, namely, 
Li doping of graphene, which has been proposed as a suitable hydrogen storage system We 
use a similar unit cell as for the boron sheet, containing six carbon atoms and three hexagons, one 
of which contains a Li atom. After relaxation, the lattice parameter of the graphene— Li cell is 4.28 
A, whereas the boron- Li cell parameter is slightly larger (5.08 A). We find that the binding of Li 
atoms to graphene is considerably weaker than the binding to the boron sheet. On graphene, the 
binding energy E\> = —0.61 eV. Obviously, a negative value of the binding energy indicates that 
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(a) 





1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 

Reaction Coordinate 



Figure 4: Diffusion of a Li atom over a (a) boron sheet and (b) graphene. The lowest energy 
diffusion paths are shown in top and perspective views. Calculated energy profiles along the paths 
are plotted on the right-hand side. The yellow dots in the plot correspond to the positions shown 
on the left-hand side. The diffusion barriers are 0.45 and 0.25 eV on the boron and graphene, 
respectively. 

this system is thermodynamically unstable with respect to bulk Li and graphene, in contrast to the 



boron-Li system, cf. Table 1 



Moreover, a Bader analysis carried out on the graphene-Li system reveals that only 0.64e is re- 
moved from Li atoms after deposition, which is significantly smaller than the case of the boron— Li 
system. These results suggest that the boron sheet is a more promising candidate for lightweight 
AM doping than pure graphene. 



Mobility of Alkali Metal Atoms 

The stability of a material such as a doped boron or graphene sheet also depends on the mobility of 
the AM atoms on the sheet. Within the transition state theory model, the most important parameter 
determining the mobility is the energy barrier for diffusion. We determine the diffusion path of a 
single Li atom on the boron sheet using the nudged elastic band method.^ 1 I 52 I 53 I The calculations 
are carried out using a 2 x 2 supercell with 32 boron atoms, where one empty hexagon is occupied 
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by a Li atom. A string of 17 images is chosen to model the diffusion path between two empty 



hexagons, as shown in Figure 4i. Positions over the centers of the empty hexagons represent 
absolute energy minima, and positions over the boron atoms at the centers of the other hexagons 
represent secondary adsorption minima. An optimal path runs over B atoms, and along B— B 



bonds, as shown in |Figure 4ft , with maxima over B— B bond centers. The calculated diffusion 
barrier is 0.45 eV. The lowest energy path for diffusion between two energy minima corresponding 
to empty hexagons proceeds via a secondary minimum over a filled hexagon, giving rise to a curved 
path. The simpler direct path between two empty hexagons involves an energy barrier of 0.55 eV 
and is therefore less likely. 

We can compare this to diffusion of a Li atom over graphene. Positions over the centers of 
the carbon hexagons represent the absolute energy minima. A Li atom diffuses over a C— C bond 



following a simple path, with a maximum over a C— C bond center, as shown in Figure 4 ~>. The 
calculated diffusion barrier is = 0.25 eV. Consequently, Li atoms on graphene are quite mo- 
bile. Using an Arrhenius expression for the jump rate v = VoexplE^/kT] with a typical attempt 
frequency of Vo ~ 10 12 Hz, one obtains v ~ 10 8 Hz at room temperature and v ~ 1 Hz at T = 105 
K. Because of their negative binding energy with respect to bulk Li, as discussed in the previous 
section, there is a thermodynamic driving force for Li atoms to cluster, which would seriously 
hamper the use of doped graphene as a hydrogen storage material.^^ It is, however, possible to 
use graphene as a template for polylithiated molecules to create a stable storage material.^ 

The diffusion barrier, E$ = 0.45 eV, for Li atoms on the boron sheet is substantially larger 
than that on graphene. Using the same Arrhenius expression as above, a jump rate of v ~ 1 Hz is 
obtained at T = 190 K. At room temperature, individual Li atoms are still mobile (v ~ 10 4 Hz), but 
the formation bulk Li metal is thermodynamic ally unfavorable because of the high binding energy 
of the Li atoms on the boron sheet. The energy difference between the minimum energy sites of a 
Li atom (over the centers of empty boron hexagons) and the secondary minima is 0.3 eV. It means 
that at room temperature only the minimum energy sites are occupied. Since the Li atoms carry a 
substantial effective charge, there is in addition no driving force for clustering of metal atoms on 
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Table 3: Calculated Binding Energies of Alkali Metals on a Boron Sheet Where Both Sides of the 
Sheet are Functionalized ( |Figure 3 ), with Respect to Bulk Metal (2?b) and Isolated Atoms (E at ), 



Equilibrium Distances to the Boron Plane (d), Bader Charges on the Alkali Atoms (0, Hydrogen 
Binding Energies (En 2 ), an d Gravimetric Hydrogen densities of the BsAM2 Systems. 



AM £ b (eV) 


£ a t(eV) 


d(A) 


Q(e) 


£ H2 (eV) 


H 2 (wt %) 


Li 0.18 


1.80 


1.62 


+0.82 


0.35-0.15 


10.75 


Na 0.02 


1.11 


2.28 


+0.60 


0.13-0.07 


8.36 


K 0.35 


1.22 


2.75 


+0.47 


0.06 


2.39 



the sheet. The minimum energy sites on the boron sheet form a regular lattice with a spacing of 
5.08 A between the lattice points, which forms the ideal template for creating a stable lattice of Li 
atoms. 

Hydrogen Adsorption on Alkali Doped Boron Sheets 

Starting from the AM lattice on the boron sheet, as discussed in the previous sections, we find in 



fact that both faces of the boron lattice can be functionalized with AMs, as depicted in Figure 5 
The chemical formula of the AM doped boron unit cells is then BsAM2 and the optimized lattice 
constants are 5.10, 5.08, and 5.08 A for the AM = Li, Na, and K, respectively. The optimized 
structures have mirror symmetry with respect to the boron plane. The binding energies of the 



AM atoms, their distances to the boron plane, and their Bader charges are given in Table 3 A 



comparison to the data of Table 1 demonstrates that there is very little interaction between the AM 
atoms located above and below the boron sheet. Our general observations for AM doping remain 
valid; i.e., also for the doubly functionalized boron sheet, we have a stable lattice of effectively 
charged AM atoms. 

Electrostatic and polarization interactions contribute most to the binding between the AM 
atoms and hydrogen molecules.^ One may therefore expect that the effective charges on the metal 
atoms after deposition onto the boron sheet determines the strength of the interaction between the 
host material and the hydrogen molecules. These charges decrease along the series Li, Na, and K, 



as shown by Table 3 The binding energies of hydrogen molecules to the AM doped boron sheet 



follow this trend. The first two H 2 molecules bind to BgLi2 with a binding energy of 0.35 eV/H2 
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(b) 





Figure 5: Optimized structure of a boron sheet doped with Li atoms on both sides, shown from (a) 
the top and (b) the side. Boron sheets doped with other lightweight AMs have a similar geometry. 




Figure 6: The optimized geometry of the Li doped boron sheet after full hydrogenation, shown 
from (a) the top and (b) the side. Three H2 molecules surround each Li atom. 



(one H2 binds to the top Li, and one to the bottom Li). The corresponding numbers for Na and K 
are significantly smaller, i.e., 0.13 and 0.06 eV/H2, respectively. 

Each Li atom in this structure can capture up to three H2 molecules, but the average binding 
energy then drops to 0.15 eV/EL-. In the fully loaded case, the BgLi2 compound then contains 10.7 
wt % hydrogen, and the chemical formula of the unit cell is BsLi2(H2)6- The optimized geometry 
is shown in |Figure 6 The centers of all the H2 molecules surrounding the Li atoms are located at an 
equivalent distance of 2.35 A from the metal atoms. These results are in line with the interactions 
of H2 molecules with isolated AM cations.ESEl] The maximum number of H2 molecules that can 
be gathered around a metal center is determined by the Coulomb interactions that develop between 
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Figure 7: Electron localization function of the boron— Li system at its fully hydrogenated state. 
The electrons are localized at the H2 molecules and around the boron sheet, where the bonding is 
covalent. The positively charged Li ions are depleted from electrons. 



Li— H2 and H2— H2, and by steric effects. 

The BgNa2 system can also capture three H2 molecules per metal. The average binding energy 
for the H2 molecules then drops to 0.07 eV/H2, which is still larger than the binding between a 
H2 molecule and a bare boron sheet. In contrast, the B8K2 system is capable of capturing only 
a single H2 molecule per K atom with a low binding energy of 0.06 eV/H2. It is clear that, if 
compared to the Li doped system, the higher electron densities on the Na and K atoms result in 
weaker interactions with H2 molecules. The Li doped boron sheet is the best option for hydrogen 
storage. 

It is interesting to note that the amount of charge on the Li atoms, as a consequence of their 
interaction with the boron lattice, stays almost constant during the consecutive addition of H2 
molecules. The H— H bonds of the adsorbed H2 molecules are only elongated by ~ 1.5 % as 
compared to an isolated H2 molecule, which illustrates that the interaction between H2 molecules 
and the Li atoms is mostly electrostatic and not chemical bonding. 



The calculated electron localization function for the BsLi2(H2)6 system is shown in Figure 7 



The Li atoms are depleted, consistent with our Bader analysis. Hydrogens stay in molecular form 
and are surrounded by their localized electron clouds. There is no direct indication of an orbital 
interaction between the Li atoms and the hydrogens. Moreover, the pattern of the covalent bonding 
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within the boron lattice is not significantly altered by the Li atoms or the H2 molecules. 

Conclusions 

We study the planar boron a-sheet as a physisorption template for hydrogen storage. Direct ph- 
ysisorption of H2 on a boron sheet gives a binding energy of 0.05 eV. Alkali metal (AM) atoms 
can be dispersed on the boron sheet to increase the hydrogen binding energies. The boron sheet 
interacts strongly with AM atoms, and the doped systems are thermodynamically stable with re- 
spect to clustering of metal atoms. Moreover, the unique geometry of the boron plane provides a 
natural lattice for the metal atoms whose nearest neighbor distance can be tuned by varying their 
concentration. In particular, Li is found to be a promising doping element for hydrogen storage. 
The strong interaction between the boron sheet and the Li atoms results in an almost complete 
transfer of the Li valence electrons to the boron sheet. Electrostatic interactions between the well- 
exposed Li atoms and the H2 molecules then lead to an average binding energy of 0.15 eV/H2. 
The system physisorbs up to a maximum of 10.7 wt % hydrogen. In contrast, AM atoms dispersed 
on graphene are generally unstable with respect to the bulk metal, which suggests that AM doped 
nanostructures based upon boron are a much better starting point for hydrogen storage than their 
carbon based counterparts. 
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